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Abstract. We review some recent results for elastic electron deuteron scattering (deuteron form 
Q , factors) and photodisintegration of the deuteron, with emphasis on the recent high energy data from 

■ Jefferson Laboratory (JLab). 
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DEUTERON WAVE FUNCTIONS AND FORM FACTORS 



Calculations of deuteron form factors and photo and electrodisintegration to the A^A'^ 
O ■ final state require a deuteron wave function. The best nonrelativistic wave functions 

are calculated from the Schrodinger equation using a potential adjusted to fit the A^A'^ 
scattering data for lab energies from to 350 MeV. The quality of realistic potentials 
O I have improved steadily, and now the best potentials give fits to the NN database with a 

^ ■ %^/d.o.f ^ 1 . The Paris potential [2] was among the first potentials to be determined from 

such realistic fits, and it has since been replaced by the Argonne VI 8 potential (denoted 
Q \ by AVI 8) [3], the Nijmegen potentials [4], and most recently by the CD Bonn potentials 

^ ■ [5,6]. The momentum space S and D state wave functions determined from three of these 

models and two relativistic models (Model IIB [7] and Model W16, one of a family of 
models with varying amounts of off-shell sigma coupling introduced in connection with 
^ ■ relativistic calculations of the triton binding energy described in Ref. [8]) are shown in 

Fig. 1. The figure shows that the S and D-state components of all of these models are 
almost identical (i.e. variations of less than 10%) for momenta below about 400 MeV, 
and vary by less than a factor of 2 as the momenta reaches 1 GeV 

Elastic electron-deuteron scattering is described in the one-photon exchange approxi- 
mation by three deuteron form factors [9, 10, 1 1] . In its most general form, the relativistic 
deuteron current can be written [9, 12] 

+GM(e2)[^^(r-^)-r"(^-^)], (1) 

where ^, J (^', J') are the polarization and momentum vectors of the incoming (outgoing) 



' This talk is a shorter version of a review being prepared for Journal of Physics G [I]. 
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FIGURE 1. Momentum space wave functions for five models mentioned in the text: AVI 8 (solid), 
Paris (long dashed), CD Bonn (short dashed), IIB (short dot-dashed), and W16 (long dot-dashed) The 
wave functions in the right panel have been divided by scaling functions for easy comparison (see Ref. 
[1] for details). 



deuterons, and the form factors G,(2^), / = 1 — 3, are all functions of = —q^, the 
square of the/owr-momentum transferred by the electron, with q = d' — d.ln practice, 
Gi and G3 are replaced by a more physical choice of form factors 

Gc = Gi + ^r[GQ 

Ge = Gi-GM + (l+ri)G3, (2) 

with T] = Q^/4m^. At = 0, the form factors Gc, Gm, and Gq give the charge, magnetic 
and quadrupole moments of the deuteron 

Gc(0) = 1 (in units of e) 

Gq{0) = Qd (in units of e/m^) 

Gm(0) = i^d (in units ofe/2m^) . (3) 

The structure functions A and B, and the polarization transfer coefficient T20 depend 
on the three electromagnetic form factors 

A{Q') = GliQ') + l^'GliQ^) + l^GliQ') 

B{Q^) = ^ni\+^)Gl{Q^) 

,4) 

where y = 2riGg/3 Gc, and 72o is T20 with the magnetic contributions removed. 
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FIGURE 2. The structure function A for five nonrelativistic models using the MMD nucleon form 
factors. The models are labeled as in Fig. 1. The right panel shows data and models divided by a "fit" 
described in Ref . [ 1 ] . The data are fully referenced in [ 1 ] . 



Comparison of nonrelativistic theory to data 

In the nonrelativistic theory, without exchange currents or (v/c)^ corrections, the 
deuteron form factors are 



Gc = 


G'eDc 


Gq = 


G'eDq 


Gm = 




G, 


2mp 





^mDm + G^^De 



(5) 



where G;^ and G;^ are the nucleon isoscalar form factors, and the Ds are the body form 
factors. All are functions of Q^. Hence the study of deuteron form factors is complicated 
by the fact that they are a product of the nucleon isoscalar form factors and body form 
factors. The dependence of the deuteron form factors on older models of the nucleon 
form factors is well discussed in Ref. [13]. A year ago the model of Mergell, Meissner 
and Drechsel [14] (referred to as MMD) gave a good fit, and could have been adopted 
as a standard (the new Gep/Gmp data from JLab [15] may change this view). The 
high data for A provide the most stringent test. In Fig. 2 we compare the data for 
A with nonrelativistic calculations using the five nonrelativistic wave functions shown 
in Fig. 1. The calculations use Eq. (5) with MMD isoscalar nucleon form factors and 
nonrelativistic body form factors. In the right panel the data and models have been 
divided by the "fit" described in Ref. [1]. 

It is easy to see that the nonrelativistic models are a factor 4 to 8 smaller than the 
data for > 2 GeV^. Furthermore, since the difference between different deuteron 
models is substantially smaller than this discrepancy, it is unlikely that any realistic 
nonrelativistic model can be found that will agree with the data. If the nucleon isoscalar 
charge form factor were larger than the MMD model by a factor of 2 to 3 it might 
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FIGURE 3. Exchange currents that might play a role in meson theories, (a) Large I — 1 7t,p, and 
A currents that do not contribute to the deuteron form factors, and (b) possible 7 = currents that are 
identically zero. The currents that do contribute to the deuteron form factors are shown in the second row: 
(c) "pair" currents from nucleon Z-graphs; (d) "recoil" corrections; (e) two pion exchange (TPE) currents; 
and (f) the famous pjty exchange current. 



explain the data, but this is also unlikely since the variation between nucleon form factor 
models is substantially smaller than this. We are forced to conclude that these high 
measurements cannot be explained by nonrelativistic physics and present very strong 
evidence for the presence of interaction currents, relativistic effects, or possibly new 
physics. 



Relativistic calculations and new physics 

The differences between the data and the nonrelativistic theory can only be explained 
by a combination of the following effects 

• interaction (or meson exchange) currents; 

• relativistic effects; or 

• new (quark) physics. 

The only possibilities excluded from this list are variations in models of the nucleon 
form factors, or model dependence of the deuteron wave functions. Previously we have 
argued that neither the current uncertainty in our knowledge of the nucleon form factors, 
nor the model dependence of the nonrelativistic deuteron wave functions is sufficient to 
provide an explanation for the discrepancies. 

Possible interaction currents that might account for the discrepancy are shown in 
Fig. 3. Because the deuteron is an isoscalar system, the familiar large 7=1 exchange 
currents are "filtered" out and only 7 = exchange currents can contribute to the form 
factors. The 7 = currents tend to be smaller and of a more subtle origin. The nucleon 
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FIGURE 4. The relativistic decision tree discussed in the text. 
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Z-graphs, Fig. 3c, and the recoil corrections, Fig. 3d, are both of relativistic origin. 
(The recoil graphs will give a large, incorrect answer unless they are renormalized 
[16, 17, 18].) The two-meson exchange currents should be omitted unless the force also 
contains TPE forces. The famous pTiy exchange current is very sensitive to the choice 
of pny form factor, which is hard to estimate and could easily be a placeholder for new 
physics arising from quark degrees of freedom. 

In most calculations based on meson theory, the two pion exchange (TPE) forces and 
currents are excluded, and, except for the pny current (which we will regard as new 
physics), the exchange currents are of relativistic origin. Additional relativistic effects 
arise from boosts of the wave functions, the currents, and the potentials, which can be 
calculated in closed form or expanded in powers of (v/c)^, depended on the method 
used. At low calculations may be done using effective field theories in which a 
small parameter is identified, and the most general (i.e. exact) theory is expanded in 
a power series in this small parameter. In these calculations, an expansion of relativistic 
effects in a power series in (v/c)^ is automatically included. Hence, with the exception of 
effective field theories, any improvement on nonrelativistic theory using nucleon degrees 
of freedom leads us to relativistic theory. 

Alternatively, one may seek to explain the discrepancy using quark degrees of freedom 
(new physics). When two nucleons overlap, their quarks can intermingle, leading to 
the creation of new A^A'^ channels with different quantum numbers (states with nucleon 
isobars, or even, perhaps, so-called "hidden color" states). These models require that 
assumptions be made about the behavior of QCD in the nonperturbative domain, and are 
difficult to construct, motivate, and constrain. At very high momentum transfers it may 
be possible to estimate the interactions using perturbative QCD (pQCD). Very little has 
been done using other approaches firmly based in QCD, such as lattice gauge theory or 
Skrymions. 

In deciding which relativistic method to use, it is first necessary to decide whether or 
not to allow antiparticle, or negative energy nucleons to propagate as part of the virtual 



intennediate state. Since nucleons are heavy and composite, so that their antiparticle 
states are very far from the region of interest, some physicists believe that intermediate 
states should be built only from positive energy nucleons, and that all negative energy 
effects (if any) should be included in the interaction. These methods are referred to col- 
lectively as hamiltonian dynamics and are represented by the left hand branch shown in 
Fig. 4. Unfortunately, it turns out that this choice precludes the possibility of retaining 
the properties of locality and manifest co variance enjoyed by field theory. Alternatively, 
in order to keep the locality and manifest covariance of the original field theory, other 
physicists are willing to allow negative energy states into the propagators. These meth- 
ods, represented by the right-hand branch of the figure, are referred to collectively as 
propagator dynamics. However, including negative energy states tends to make calcula- 
tions technically more difficult and harder to interpret physically, and those who advo- 
cate the use of hamiltonian dynamics do not believe the advantages of exact covariance 
justify the work it requires. 

Unfortunately, these two methods are so fundamentally different that many physicists 
do not realize that the limitations of one may not apply to the other. For example, for 
some choices of propagator dynamics all 10 of the generators of the Poincare group 
will depend only on the kinematics, and the Poincare transformations of all amplitudes 
can be done exactly. With hamiltonian dynamics this is not the case; some of the 10 
generators must depend on the interaction, and transformation of matrix elements under 
these "dynamical" transformations must be calculated. Comparison of the two methods 
is therefore very difficult; the language and issues of each are very different and one can 
be easily misled by the different appearance of the results. 



Comparison of relativistic calculations with data 

The high predictions for 7 relativistic models using A^A^ degrees of freedom and 
one quark cluster model are shown in Figs. 5-6. They include 

• two propagator calculations: VGO [19] (using the Spectator equation), and PWM 
[20] (using the modified Mandelzweig- Wallace equation); 

• two instant-form calculations: FSR [21] (without a v/c expansion) and ARW [22] 
(using a v/c expansion); 

• two front-form calculations: CK [23] (with the light front retained as an unphysical 
degree of freedom) and LPS [24] (using a specially constructed current operator); 

• a point-form calcualtion: AKP [25]; and 

• a quark model calculation: DB [26]. 

The model dependence of the eight calculations is large. Figure 5 shows the predic- 
tions for A (2^). In these figures we have intentionally left out the model dependent pTiy 
exchange current from all of the calculations. All of the models except the AKP point- 
form calculation give a reasonable description of A out to 2^ ~ 3 GeV^, beyond which 
they begin to depart strongly from each other and the data. Taking into account that the 
pjiy exchange current could be added to any of these models, and that this contribution 
tends to increase A above 2^ ~ 3 GeV^, four models seem to have the right general 




FIGURE 5. The structure function A for the eight models discussed in the text. Left panels show the 
propagator and instant-form results: FSR (solid line), VOG in RIA approximation (long dashed line), 
ARW (short dashed line), and PWM (dotted line). Right panels show the front-form CK (long dashed line) 
and LPS (short dashed line), the point-form AKP (medium dashed line) and the quark model calculation 
DB (solid line). In every case the calculations have been divided by a scaling function given in Ref. [1]. 
The data are labeled as in Fig 2. 



behavior: the VOG, FSR, ARW and the quark model of DB (but there are no results 
for this model beyond = 4 GeV^). None of these models fit the data without a pny 
exchange current, and all models would be improved by adding such a current contribu- 
tion, showing that there is some evidence for new physics at high Q^. Ironically, none of 
the models favored by the high data does as well at low as the three "unfavored" 
models shown in the right panels (unless the Platchkov [27] data are systematically too 
low). 

Finally, Fig. 6 shows the predictions for the structure functions A, 5, and T2Q for the 
eight models. The LPS calculation shows a large discrepancy with the T2Q data, but the 
most striking feature of these plots is the large model dependence of the predictions for 
B{Q^)- The magnetic structure function provides the most stringent test of the models. 
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FIGURE 6. The structure functions A, B, and T20 for the eight models discussed in the text. VOG full 
calculation (CIA plus pny- solid line); VOG in RIA (long dashed line); FSR (medium dashed line); ARW 
(short dashed line); DB (widely spaced dotted line); CK (long dot-dashed line); AKP (short dot-dashed 
line); PWM (dashed double-dotted line), and LPS (dotted line). The data are as labeled in Fig. 2. 




and the predictions are comparatively free of the pTry exchange current (which gives only 
a small contribution to B). Examination of the figure shows that the B predictions of the 
PWM, ARW, AKP, CK models fare the worst. In all, taking the predictions for the three 
structure functions together, the best results are obtained with the FSR, VOG, and DB 
models. 



Conclusions (deuteron form factors) 



What have we learned from our measurements of the deuteron form factors? Our 
comparison of theory and experiment leads to the following conclusions: 

• Nonrelativistic quantum mechanics (without exchange currents or relativistic ef- 
fects) is ruled out by the A(2^) data at high Q^. Reasonable variations in nucleon 
form factors or uncertainties in the nonrelativistic wave functions cannot remove 
the discrepancies. 

• In approaches using ATAT degrees of freedom only, relativistic effects (or / = 
meson exchange currents) could be large enough to explain the data. 

• Some models that include relativistic effects (or meson exchange currents) and use 
NN degrees of freedom with realistic forces are close to the data. None are entirely 
satisfactory. 

• The model dependence of relativistic effects (or meson exchange currents) is larger 
than the errors in the data, even at low Q^, and is not understood. 

• There is evidence that new physics (either in the from of the pTry exchange current 
or something else) is beginning to show up in the A structure function above of 
2 - 3 GeV2. 

• The deuteron form factors provide no evidence for the onset of pQCD, but quark 
cluster models could explain the data. 

Study of the experimental situation leads to the following conclusions: 

• The minimum of B is very sensitive to details of the models, and improved mea- 
surements of B for in the region 1.5-4 GeV^ are particularly compelling. It is 

important to accurately map out the zero in the B structure function. 

• Detailed disagreements between theories and different data sets suggests the need 
for precision studies at low 

THRESHOLD ELECTRODISINTEGRATION 

Threshold deuteron electrodisintegration measures the <i(e, e')pn reaction in kinematics 
in which the proton and neutron, rather than remaining bound, are unbound with a few 
MeV of relative kinetic energy in their center of mass system. If the final state energy 
is low enough, the final state will be dominated by transitions to the final state, and 
will be a pure AS = 1, A/ = 1, MX transition, similar to the — > A(1232) transition. 
This transition is a companion to the B structure function; both are magnetic transitions 
and both are filters for exchange currents with only one isospin (J ^ J is A/ = and 
d — > is A/ = 1). To see the similarity, compare the top right panel of Fig. 6 with the 
threshold measurements shown in left panel of Fig. 7. Both have a similar shape, and in 
both cases the uncertainties in the theoretical predictions are large. 

The similarity of these two processes (elastic and threshold inelastic) also holds for 
the theory. These two processes can be used to separately determine the precise details 
of the / = and 7=1 exchange currents. Once the exchange currents are fixed, they can 




FIGURE 7. Left panel: The cross section for threshold electrodisintegration of the deuteron. (See 
Ref. [1] for discussion of the data and the theory.) Right panel: The variation of with the photon 
energy v for various values of x. The shaded regions show the approximate thresholds for the production 
of bands of nucleon resonances. The numbers in the small circles are the number of distinct channels in 
each band. 



be used to predict the results of d{e, e'p)n over a wide kinematic region. Any theoretical 
approach that works for the form factors should also work equally well for threshold 
electrodisintegration, yet very few of the groups who have calculated form factors have 
also calculated the threshold process. These calculations, when completed, will provide 
a more definitive test of the various relativistic approaches discussed in the previous 
sections. 



DEUTERON PHOTODISINTEGRATION 

Since both the recent deuteron form factor measurements and the recent high energy 
deuteron photodisintegration measurements have been made with 4 GeV electron beams, 
it is sometimes assumed that the same theory should work for both. This need not be 
the case, because the kinematics of elastic electron-deuteron scattering and deuteron 
photodisintegration are very different, and the physics being explored by these two 
measurements is also very different. The implications of this remarkable feature of 
electronuclear physics is often not fully appreciated. 

The kinematics of elastic scattering and photodisintegration are compared in the right 
panel of Fig. 7, which shows — as a function of the photon (real or virtual) 
energy. Here x = 2^/(2mv) is the familiar Bjorken scaling variable. The mass of the 
final excited state increases rapidly as x decreases below its maximum allowed value 
of X = m^/m ~ 2. For any energy v or any Q^, elastic ed scattering leaves the pn 



system bound, with no internal excitation energy added to the two nucleons. As x — > 
(the real photon limit) the maximum value of W is reached for any given beam 
energy. The 24 well established nucleon resonances, and the bands of thresholds at 
which these resonances are excited, are show in Fig. 7. At E'y = 4 GeV, the final state 
mass is approximately 4.5 GeV, and at least 286 thresholds for the production of pairs 
of baryon resonances have been crossed (and there are probably more from unseen or 
weakly established resonances). A photon energy of 4 GeV corresponds to np scattering 
with an np laboratory kinetic energy of about 8 GeV! 

It is clearly very difficult (if not impossible) to construct a theory of high energy pho- 
toproduction in which all of these resonances and their corresponding 286 production 
thresholds are treated microscopically. By contrast, elastic electron deuteron scattering 
requires a microscopic treatment of only one channel. All of the 286 channels also con- 
tribute to elastic scattering, of course, but in this case they are not explicitly excited, and 
can probably be well described by slowly varying short-range terms included in a meson 
exchange (or potential) model. In photodisintegration, each of these channels is excited 
explicitly and an alternate framework that averages over the ejfects of many hadronic 
states is needed. The alternatives are to use a Glauber-like approach, or to borrow from 
our knowledge of DIS and build models that rely on the underlying quark degrees of 
freedom. 



High Energy Photodisintegration 

Figure 8 shows the published high energy photodisintegration data, from experiments 
NE8 [28, 29] and NE17 [30] at SLAC, and E89-012 [31] and E96-003 [32] at CEBAF. 
These experiments determine cross sections for 9cm ~ 36°, 52°, 69°,and 89° at energies 
from about 0.7 to 5.5 GeV; there are also some backward angle data up to 1.6 GeV from 
NE8. 

The main feature of the data above about 1 GeV is the s~^^ (s~^^) fall off (where 
s = {pi+ PiY is the square of the cm. energy) of the cross sections da/dt (do/dQ.) 
at Gem = 90 and 69°, in agreement with perturbative QCD expectations. In contrast, the 
cross sections at the forward angles 36 and 52°, fall off more slowly, with s"^ scaling 
at lower energies, until the onset of the 5 ^ ^ behavior at about 4 and 3 GeV beam energy, 
respectively. At each angle the onset of the s~^^ behavior corresponds to a perpendicular 
momentum, pr, of approximately 1 GeV^. 

The highest energy polarization measurements are of py (the induced polarization of 
the proton), and Cy and Q' (the transfer of circular polarization from the photon to the 
proton) from CEBAF E89-019 [33]. The left panel of Fig. 9 shows the striking feature 
that the induced polarization py is consistent with zero (as predicted by pQCD) at ener- 
gies above about 1 GeV, the same energy at which the s~^^ cross section scaling begins. 
The right panel of the same figure shows that the polarization transfer observables both 
appear to peak near 1 GeV, and decrease at higher energies. 
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FIGURE 8. Photodisintegration cross section s^^da/dt versus incident lab photon energy. The calcula- 
tions are from Kang, Erbs, Pfeil and Rollnik (solid line), Lee (dashed line), Raydushkin quark exchange 
(dot-dashed line), reduced nuclear amplitudes of Brodsky Hiller (dotted line), quark qluon string model 
(short dashed line), and Frankfurt, Miller, Strikman and Sargsian QCD rescattering (shaded region). (See 
Ref. [1] for references and discussion of the data and theory.) 



Conclusions (deuteron photodisintegration) 

Review of deuteron photodisintegration suggests the following: 

• A microscopic meson-baryon theory of deuteron photodisintegration must describe 
the A^A^ interaction at high energies, including pion production and the contributions 
of hundreds of A'^* channels. It is unlikely that such a theory will be constructed in 
the foreseeable future. 
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FIGURE 9. Left panel: Induced polarization for deuteron photodisintegration at 9cm = 90°. The cal- 
culations are from Kang, Erbs, Pfeil and Rollnik (dashed line) and from Sargsian (dot-dashed line). 
Right panel: Polarization transfer for deuteron photodisintegration at 9cm = 90°. The calculation is from 
Schwamb and Arenhovel. 



• For pj>l GeV^, cross sections appear to follow the constituent counting rules, 
but it is expected that an absolute pQCD calculation would greatly underpredict the 
data. Similar observations may be made for other photoreactions, and it remains to 
be seen how this behavior arises, and if there is a general explanation for it. 

• Some nonperturbative quark models do well describing the data qualitatively. Fur- 
ther theoretical development and experimental tests of these models would be de- 
sirable. 



OVERALL CONCLUSIONS 

Our overall conclusions from the study of form factors and high energy photodisintegra- 
tion can be briefly summarized as follows: 

• Meson theory works well in cases where all of the active hadronic channels that 
can contribute to a process are included. This has been done for the deuteron form 
factors (where only the A^A^ channel is active), but is impossible for high energy 
deuteron photodisintegration where lOO's of N*N* channels are active. At high 
energy, any successful meson theory must include relativistic effects. 

• New approaches, probably using quark degrees of freedom, are needed for high 
energy deuteron photodisintegration. 

• Meson theory behaves as might be expected. It works for the deuteron form factors 
and does not work for photodisintegration (because of the number of channels). 

• Perturbative QCD does not behave as might be expected. Theoretically it should 
work (or not work) equally well for the deuteron form factors at high as it 
does for photodisintegration at high (both are exclusive processes). However 



the scaling laws seem to work qualitatively for photodisintegration and to fail badly 
for elastic scattering. 

Please refer to Ref. [1] for a more complete discussion of all of the points in this talk. 
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